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We investigate the phase structure of massless three-flavor QCD by extending the Nambu-Jona- 
Lasinio model to include the effects of confinement and the axial anomaly. We study the interplay 
between the chiral and diquark condensates induced by the axial anomaly, as well as their relation- 
ship with the Polyakov loop, which parameterizes confinement. By minimizing the thermodynamic 
potential we construct the QCD phase diagram and investigate the possibility of realizing a recently 
discovered low temperature critical point and an associated BEC-BCS crossover. We also perform 
a Ginzburg-Landau expansion of the thermodynamic potential, comparing our results to a prior 
analysis based purely on symmetry considerations, in order to assess the lowest-order effects of the 
condensate-confinement couplings. 



I. INTRODUCTION 

The phase structure of strongly-interacting matter has 
been a subject of great interest since the emergence of 
quantum chromodynamics (QCD) in the early 1970s. 
However, it has gained greater attention in recent years as 
the boundaries of experimental study have continued to 
expand, particularly at the Relativistic Heavy Ion Col- 
lider (RHIC), and even now at the Large Hadron Col- 
lider (LHC). In addition, more powerful computational 
techniques continue to produce more reliable lattice data, 
which may be compared with experimental results. 

However, there remains much wanting in the theo- 
retical analysis of QCD matter, due to the dual prob- 
lem of the intractibility of nonperturbative QCD and 
the fermion sign problem, which limits our ability to 
extend lattice techniques to non-zero chemical poten- 
tial. One tool which has proven useful in filling this gap 
is symmetry-based models, which have been developed 
to describe critical characteristics of strongly-interacting 
matter including chiral symmetry breaking, diquark pair- 
ing, and confinement, while still retaining calculational 
tractibility [H-Ql- Recently, increased attention has been 
given to the importance of the instanton-induced ax- 
ial anomaly and the attractive chiral-diquark condensate 
coupling that it mediates 

The axial anomaly, which gives rise to the Kobayashi- 
Maskawa, 't Hooft (KMT) effective six-quark interaction, 
plays a crucial role in determining both the properties 
of the light mesons and the phase structure of quark 
matter. In the first case, the anomaly is responsible 
for breaking the U(l)^ symmetry of QCD and giving 
rise to the anomalously large mass of the rj' meson. In 
the second, recent work by Hatsuda et al. has demon- 
strated the importance of the axial anomaly in determin- 
ing the topology of the QCD phase diagram I, In 
particular, under the right circumstances the anomaly- 
induced attraction between chiral and diquark conden- 
sates leads to a low temperature (T) critical point and a 
corresponding BCS-BEC crossover between the chirally 




FIG. 1: Schematic QCD phase diagram showing regions of bro- 
ken chiral symmetry (hadron), BCS quark pairing and color su- 
perconductivity (CSC), coexistence of chiral and diquark con- 
densates (COE), and a deconfined quark-gluon plasma (QGP). 
Single lines denote second order phase transitions while double 
lines denote first-order transitions. Adapted from [lo| . 



broken Nambu-Goldstone (NG) phase and a color super- 
conducting (CSC) phase characterized by diquark pair- 
ing, as shown in Fig. [T] 

One useful method for describing dense quark matter 
is the Nambu-Jona-Lasinio (NJL) model |ll| - [l7j . First 
developed to describe chiral symmetry breaking, it has 
been extended to include diqtiark pairing and confine- 
ment, via the Polyakov loop [18l - [25| . The resulting PNJL 
model has been shown to possess coincident chiral and 
deconfinement transitions for both two- and three-fiavor 
systems, in close agreement with current lattice simula- 
tions [26-29j. More recently, the NJL model has also 
been used to investigate the effects of the axial anomaly 
on the QCD phase diagram fll| . 

In the present paper we further develop recent exten- 
sions of the NJL model by including the effects of both 
the axial anomaly and confinement on the massless three- 
flavor QCD phase diagram. In Sec.|ll]we review the NJL 
model including the KMT interaction, while in Sec. IIIII 
we introduce the Polyakov loop and discuss its ability 
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to describe confinement. In Sec. IIVI we combine these 
two models and derive the thermodynamic potential of 
the three-flavor PNJL model. In Sec. |V] we construct 
the QCD phase diagram by minimizing the thermody- 
namic potential, paying special attention to the recently 
discovered low T critical point and the criteria for its 
emergence. Finally, in Sec. IVII we perform an expansion 
of the thermodynamic potential of the Ginzburg-Landau 
form and assess the lowest-order couplings between the 
quark condensates and the Polyakov loop. 



II. THREE-FLAVOR NJL MODEL WITH AXIAL 
ANOMALY 

Following Abuki et al. we write the three-flavor NJL 
Lagrangian 

Cnjl = q{i$ -m + ^if )q + £(4) + £(6)^ (i) 

where q — (u,d,s)"^ is the quark field, to = 
diag(TOi,, TOc(, TOs) is the quark mass matrix in flavor 
space, which we choose to be diagonal {rh = toI), is 
the quark chemical potential, and C^^^ and £^^) are four- 
and six-quark interaction terms respectively. 

The standard form of the four-quark interaction C^^^ 
is chosen to respect the chiral symmetry of QCD and 
contains color-flavor-locked (CFL) quark-quark interac- 
tions [1,[33] 

£(4)^44)^^(4)^ (2) 
8 

£i^^ = GY,[{qTaqf + {qij5raqf, (3) 
a=0 

+ [qrA\A'Ct){q^CTA\A-q)], (4) 

where are the U(3) flavor generators (a = 0...8), ta 
and \a are the antisymmetric flavor and SU(3) color gen- 
erators {A, A' = 2, 5, 7), and C is the charge conjugation 
operator. The generators are normalized by the condition 
Tr(TaTb) = 25ab^ and we take G, H > 0, corresponding to 
attractive interactions. We also find it useful to define 



the chiral and diquark operators 

<l>^, - {qR)i{qL)l (5) 

{dL)ai = ^abceijkiqL)iCiqL)c, (6) 

{dR)a^ = <^abceijkiqR)lC{qR)c, (7) 



where i,j,k and a,b,c are flavor and color indices re- 
spectively and qL and qn denote states of left- and right- 
handed chirality. In terms of these operators, the four- 
quark interactions can be written in the form 

44)=8GTr/(0t0), (8) 

=2i/Tr/(4di+4dfl). (9) 



The six-quark interaction can also be written as the 
sum of two terms 

/:(«) = 4^)+4^j, (10) 

= -8if(det0 + h.c.), (11) 
4d = K' (Tr,,/[(4di)0] + h.c.) . (12) 

The first term, c'^\ is the standard KMT interaction, 
which is the result of the instanton-induced QCD axial 
anomaly [1, 0] . The second term, c'f^ , is the effective in- 
teraction between the chiral and diquark fields, mediated 
by the QCD instanton 

The condensates favored by C^"^^ and C^-^^ are the 
flavor-symmetric chiral and diquark condensates in the 
spin-parity 0^ channel: 

{qWa) = cr*y {q^C^f,TA\A'q) = d5AA'- (13) 

In mean field the Lagrangian becomes 

Cnjl = q{i$ - M + ii-i^^)q ~ ^{A*q^Cj5TA\Aq + h.c.) 

-V{a,d), (14) 

where the sum over A — 2, 5, 7 is implied and where the 
effective quark mass is 

M = m-4Ga + 2Ka'^ + ^\d\^, (15) 

and the pairing gap A is 

A = -2d(^H - ^(T^ . (16) 

In addition, the condensates contribute to the system's 
potential directly an amount 

V{<J, d) = GGct^ + 3H\d\^ - AKcr^ — \d\^(j. (17) 

In order to diagonalize the Lagrangian, it is convenient 
to define the Nambu-Gor'kov field 

where q'-^ = Clf" is the charge-conjugate quark field. 
With this, we write the mean field Lagrangian in the 
form Cnjl = ^S~^^ — V, where the inverse propagator 
in Nambu-Gor'kov space is 

s-\k) = (^\f:^"-'' /^^"(f^^,,V (19) 

^ ^ y -A 75TaAa f - fJ-r - M J 

Having diagonalized the Lagrangian, we compute the 
thermodynamic potential by performing the Gaussian in- 
tegrals over the fields ^ and ^. Thus, we write the ther- 
modynamic potential in the standard way 

nNjL^V{a,d)-^Y. J (0j^trln[/35-i(^„,k)], 

(20) 
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TABLE I: Couplin g co nstants and dynamical quark mass for 
the PNJL model fTl|,TTl. 









M (MeV) 


1.926 


1.74 


12.36 


355.1 



where (3 = 1/T is the inverse temperature, the trace is 
taken over Dirac, color, flavor, and Nambu-Gor'kov in- 
dices, the factor of 1/2 accounts for the double-counting 
inherent in the Nambu-Gor'kov formalism, and aj„ ~ 
i{2n + l)7rT are the fermionic Matsubara frequencies. 

The couplings G, H, and K, as well as the cutoff A 
are fixed by fitting to experimentally determined mesonic 
properties, as discussed by Buballa [l^, and are given in 
TableHl The second anomaly coupling K' will be adjusted 
"by hand" to determine its effect on the low T critical 
point. In the following calculations we assume that the 
current masses of the quarks are zero. 



III. CONFINEMENT AND THE POLYAKOV 
LOOP 

Having described the portion of the Lagrangian re- 
sponsible for chiral symmetry breaking and diquark pair- 
ing, we next extend our model to include confinement. 
We begin by introducing a temporal gauge field a-nd 
replacing the derivatives in Eq. (|14p with the covariant 
form 



(21) 



where ^o(x,t) = Aq, a matrix in color space, is chosen 
to be a time-independent homogeneous gauge field. 

As shown by Polyakov, in the absence of quarks the 
order parameter for confinement is the thermal Wilson 
line with periodic boundary conditions (also called the 
Polyakov loop) [l^, [1^ Hlj . Unfortunately, despite the 
ubiquity of this object, its notation in the literature is 
very inconsistent. Thus, we define precisely the quantity 
we will use to describe the QCD deconfinement transi- 
tion. 

The gauge-invariant Wilson loop is defined as 



(22) 



L(x) =Vcxp<i (p dx^'A 



where V is the path-ordering operator. Note that £(x) 
inherits any non-Dirac indices of the gauge field A^. 
Thus, in the present theory, L(x) is a matrix in color 
space. When working at finite temperature, we perform 
a Wick rotation to Euclidean time by defining r = it, 
with the Aq component of the gauge field also rotated to 
A4 = iAo. In the absence of the spatial components of 
Afj_, the Polyakov-loop matrix is now 



L(x)=7'exp|i J dTAo(x, r) 



(23) 



where V is now understood to be the ordering operator 
in imaginary time. Specializing to constant Aq, Eq. ([ 
reduces to 



L = e 



(24) 



This is the expression that we will use throughout our 
analysis. Note that Aq is a Hermitian field, in contrast 
to some authors who write L = e^^°, with ^0 anti- 
Hermitian 0,|32. 

In the limit of infinitely massive quarks, the free energy 
of a static quark-quark pair separated by x is [13, 



g~/3F„(x) ^ ($(x)$(0)). 



(25) 



where $ and <& are the thermal expectation values of the 
traced Polyakov loop and its conjugate 

$(x) = i- (tr,L(x)) $(x) = i- (tr,it(x)) .(26) 

As |x| — ^ 00, the expectation values of the operators can 
be factored to obtain 



;<i>(x)) ($(0)) = |$(o) 



(27) 



For confined quarks, as the separation between two 
quarks tends to infinity, so does the energy required to 
maintain the separation {Foo = 00). Thus, $ is an order 
parameter for the deconfinement transition: 

$ = confined, 
$ 7^ deconfined. 

From the normalization adopted in Eq. (1261) we see that 
"complete deconfinement" is the limit $ — )• 1, while < 
$ < 1 describes a system somewhere between complete 
confinement and complete deconfinement. 

We next introduce a Polyakov-loop potential U{^,'^) 
which describes the deconfinement transition in the pure- 
gauge sector. Different forms have been proposed for 
U{^, $), but in all cases, the functional form is chosen to 
be consistent with the Z(3) center symmetry of QCD [Tol - 
[2H . In this study we follow Fukishima and RoCner and 
use the potential [l^, dH 

^ --a{T)^^ + b{T) ln[l - 6$$ 

+4($3 + $^)_3($$)2]^ (28) 

where the temperature-dependent coefficients have the 
form 

a(T) = ao + ar(f)+a,(|)' 6(T) = ^3 

(29) 



TABLE II: Coefficients of the Polyakov-loop potential pi 



ao 


ai 


0.2 


&3 


3.51 


-2.47 


15.2 


-1.75 
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FIG. 2: (color online) Dimensionless Polyakov loop potential, 
U{^)/T'^ , which models confinement in the pure-gauge sector. 
Below To the potential is minimized for $ = and the system 
is confined, spontaneously breaking QCD's Z(3) center symme- 
try. Above To the symmetry is restored and the system has a 
non-zoro <I>, indicating movement toward a quark- gluon plasma 
(QGP). 



the tti and bi are fixed by comparison with lattice data 
at fi = 0, and are shown in Table HIl plj . The parame- 
ter Tq, which is 270 MeV in the pure-gauge sector, will 
be fixed to reproduce the three-flavor lattice transition 
temperature |33l - l36| . Note that, as shown in Fig. [21 by 
diverging in the limit $ — >■ 1, the logarithmic form of 
(<&,¥) ensures that $ < 1. 



IV. THE PNJL MODEL : CONFINEMENT AND 
CHIRAL SYMMETRY 



Having modeled both chiral symmetry breaking and 
confinement, we now combine the results of the prior 
sections to write the full Lagrangian Cpnjl = ^njl ~ 
U{^, $), with the replacement — = d^ — iA^: 



L 



PNJL 



(30) 



Note that the coupling between the quarks and the 
Polyakov loop is determined solely by the covariant 
derivative, which effectively makes the replacement fj, — > 
H + Aq. Thus, in mean field the inverse propagator in 
Nambu-Gor'kov space becomes (generalizing Eq. 



s-\k) 



^ + i^i + Ao)-f° - M A-/5ta\a 

-A*j5ta\a - (m + ^0)7° - M 



(31) 

As noted, ^0 is a matrix in color space, and may therefore 
be expressed in the form ^0 = AnXa/2. In the Polyakov 
gauge, with ^0 diagonal, we write [U Isl] 



^0 = 03^3 + 08^8, 



(32) 



where A3 and Ag are the symmetric generators of SU(3) 
color. In analogy with Eq. (|20p . the thermodynamic 



potential is then 

^PNJL = V{a,d) +U{<j)3,<j)s) 



E 



cPk. 1 

j^-^ -tT\n[PS-^{uJn,k)]Ao^rAo 



(33) 

where we have gone to the Euclidean signature by taking 
Ao — > iAo, as is required at finite temperature (39l l40l |. 
and have expressed Z^($, $) in terms of the parameters 
03 and which are related via the transformation 

^ \^il3{<i>3+<t>8/V3) _,_ gi/3(-03+08/\/3) _^ g-2i/30s/\/3" 
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(34) 



From Eqs. (|3T|) and (1321) we find that the presence 
of Aq effectively renders the chemical potential color- 
dependent. In the Gell-Mann basis, in color space we 
can write /i + iAo = diag(/ii, /i2, /i3), where 



9s 
V3 



<P8 



(35) 



2iQ 



M3 = M 



V3 



The inverse propagator in Eq. ([31]) is a 72 x 72 matrix, 
and the determinant is therefore a 72nd-order polyno- 
mial in ujn- While prior work has investigated special 
cases of this expression, including the two-flavor PNJL 
model [m and the three-flavor NJL model [ll|, in this 
study we consider the general three-flavor PNJL model. 
As a result, the evaluation of the trln in Eq. ([^(I|) is 
much more difficult, and relies on general relations in- 
volving the determinants of 2 x 2 (Nambu-Gor'kov and 
Dirac) and 3x3 (color and flavor) block matrices. 

Exchanging the order of the eigenvalue and Matsubara 
sums, we make use of the relation 

1 



ln[/3(c.„±i?,)] =ln(l + e-^^^-; 



2 ^' 



(36) 



where AEj — Ej 



Ej^'^'^, the difference in the eigenvalue 
from the noninteracting case, comes from measuring f2 
relative to the noninteracting Dirac sea. Thus, we are 
able to write the thermodynamic potential in the form 

^PNJL = V{cr, d) 
1 



E 



^(03,08) 

(27r)3 



ln(l + , 



-PE,) + _13AE, 



(37) 

Solving the characteristic polynomial of Eq. ([311) yields 
18 distinct eigenvalues, each with multiplicity 4 (2 spin 
x 2 Nambu-Gor'kov). The first 48 eigenvalues (of which 
12 are distinct) are of the form 



E 



1-48 — 




A|2± 



(38) 
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where i ^ j, Ek = -y/k^ + M^, and the two ± are inde- the roots of the polynomials F(a;„, Ek) and F{ujn, —Ek) 
pendent. The remaining 24 eigenvalues (6 distinct) are which satisfy limi;_j.oo Ej = c», where 



F[Un, Ek) = (Wn + Ek+ lll)[(jJn - Ek - /il)(w„ + Ek + M2)('^n ^ Ek - fi2){'^n + Ek + /i3)(Wn - Ek - ^3) 



A|2(a;„ + 


Sfe + 




+ Ek 


+ M2)(Wn - 


Ek - 


fJ-sW^^n - Ek 


-Mi) + 


{(^71 


-Ek 




A|2(w„ 4- 


Sfe + 


A*l)(^n 


+ Ek 


+ Ai3)(Wn - 


Ek - 


/i2)[(w„ - Ek 


-Mi) + 




-Ek 


-M3)] 


A|2(w„ + 


Sfe + 




+ Ek 


+ Ai3)(Wn - 


Ek - 


/il)[(w„ - Ek 


-M2) + 




-Ek 


-Ma)] 


A|4[K- 


^Ek^ 


-Mi) + ( 


UJn + 


Sfc + M2) + 




-Ek+m)] 










x[(a; 




« - Ml) - 




- Ek- 


+ i^n - Ek- Ms)] ^ 


-4|A|^. 









The function F{LOn, Ek) is an even sixth-order polynomial 
and therefore reducible to a cubic polynomial, whose so- 
lutions can be obtained exactly. Thus, all of the model's 
eigenvalues can be obtained explicitly, and the thermo- 
dynamic potential computed via Eq. ([37]) • In order 
to construct the phase diagram, we must now minimize 
0(cr, d, 03, ^s) with respect to all variables. 

We must first ensure that the minimization of ft is well- 
defined. In fact, from Eqs. (|55)) and (15^]) we see that 
the eigenvalues are manifestly complex, and the imagi- 
nary part of f2 need not vanish. In general, fl is there- 
fore a complex function, the minimization of which is 
ill-defined. Indeed, it is not at all clear how one would 
interpret a complex potential in this context. However, 
as noted by Weise et al. in the context of the two-flavor 
PNJL model, if we assume that <f>^ and are real then 
it follows that (f)^ = so that only 03 is left as an inde- 
pendent variable [2l|. Then Eq. reduces to 
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FIG. 3: Phase diagram for the three-flavor NJL model (no con- 
finement). Thick lines represent first-order phase transitions 
while thin lines represent second order transitions. The dotted 
vertical line is BEC-BCS crossover defined by M{fi, T) = /x. The 
low T critical point is at {/i, T) = (291 MeV, 35 MeV). 



1 + 2 cos 



(40) 



For 08 = we also find ^[ = /i2 and /X3 = ^, so the 
eigenvalues come in conjugate pairs. This ensures that 
the thermodynamic potential is in fact real and that the 
minimization procedure is well-defined. 



V. PHASE STRUCTURE 

Having obtained the thermodynamic potential, we now 
construct the 3-flavor PNJL model phase diagram by 
minimizing Eq. p7l) with respect to cr, d, and 03 at each 
point (/i, T). In order to assess the effects of confinement 
on the QCD phase structure and the low T critical point 
we first construct the phase diagram for $ = 1 (no con- 
finement), and then compare to the results for the full 
PNJL model. 



A. Without Confinement 

In the absence of the Polyakov loop = 03 = 0, so 
that we set $ = 1 (see Eq. pp]) ) and eliminate W($,<I>) 



from the thermodynamic potential. The thermodynamic 
potential then reduces to 

3 



NJL 



iln(l + e-^^?) 



- ln(l + e-l^^^) + A(3/\Ef + i ^A^;, 



where the eigenvalues are now 



Ef 
Et 



v/(i?fe±M)' + |A|2 
v/(i?fe±M)^+4|A|2, 



(41) 

(42) 
(43) 



and where ^_|_ denotes summation over E+^ and E^^. 

The minimization of £7 with respect to a and d yields 
the phase diagram shown in Fig. |31 The critical tem- 
perature at /.t = is found to be T^'"" ~ 153 MeV, in 
agreement with Abuki et al., and within the margin of 
error of current lattice results [ill IHj • Note that this 
temperature is determined entirely by the couplings G, 
H, and K , which were in turn fixed by emperical mesonic 
properties {T^"^^ proves independent of K' , which de- 
scribes condensate coupling, as the diquark condensate 
is absent in this portion of the phase diagram). 
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FIG. 4: (color online) Contour plot for the chiral condensate in 
the NJL model, showing contours of c = 0.2, 0.4, 0.6, 0.8 and 1.0 
CO, where tro is the maximum value of the chiral condensate. 
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FIG. 5: (color online) Contour plot for the diquark condensate 
in the NJL model, showing contours of d = 0.2, 0.4, 0.6, 0.8 and 
1.0 do, where do is the maximum value of the diquark conden- 
sate. 



As reported by Abuki et al, the topology of the phase 
diagram depends critically on the ratio of anomaly cou- 
plings K — K' /K [ll|. We find, in agreement with their 
report, that for k < 4.2, the transition out of the Nambu- 
Golstone (NG) phase is first-order for all temperatures, 
while for k > 4.2, the low-T critical point emerges, below 
which there is a smooth crossover from the NG phase to 
a CSC-like COE phase. 

In a result not previously reported, we find that as k 
increases, the critical point moves up the NG-COE phase 
boundary (to higher T and lower /i) until it vanishes into 
the NG- "normal" (NOR) phase boundary at k = 4.8. We 
therefore find that there are three distinct structures of 
the NJL phase diagram, determined by the value of n: 
(1) K < 4.2, (2) 4.2 < K < 4.7, and (3) k > 4.8. Since we 
focus on the existence of the low T critical point, we are 
interested primarily in structure (2). Therefore, in dis- 
playing our results, we choose k — 4.2 as a representative 
value for which the critical point is realized. 

Considering the contour plot of a, shown in Fig. 01 we 
see that the chiral condensate is relatively slowly varying 
within the NG phase, which corresponds to a relatively 
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FIG. 6: Phase diagram for the three-flavor PNJL model. Thick 
lines represent first-order transitions, while thin lines represent 
second order transitions. The dotted vertical line is the BEC- 
BCS crossover, defined by M{fi, T) = fi. The low T critical point 
is at {ti,T) = (288 MeV, 36 MeV). 



constant effective quark mass. Near the phase bound- 
aries, however, a varies rapidly, either falling discontin- 
uously (NG-NOR or NG-COE, above the critical point) 
or undergoing a very rapid, though smooth, crossover 
(NG-COE, below the critical point). 

Similarly, the contours of d show that while d is not 
maximized over as great a portion of the phase diagram 
as CT, it is relatively constant throughout the COE region, 
dropping rapidly only near the second order COE-NOR 
phase transition and the NG-COE crossover (Fig. [5]). 



B. With Confinement 

Before constructing the three-flavor QCD phase dia- 
gram with confinement, as noted in Sec. IIIH we fix the 
parameter To by matching the critical temperature at 
fi — with current lattice data. With massive quarks, 
the definition of Tc in this context would be, as noted 
by Aoki et al, ambiguous, there being at least three 
standard choices: (1) a maximum of the chiral suscep- 
tibility, (2) a maximum of the quark number suscepti- 
bility, and (3) a maximum in d^/dT [35M37|. While 
these transitions are coincident in the NJL and PNJL 
models, current lattice calculations with non-zero current 
quark masses yield slightly different values for the three 
critical temperatures. However, for massless quarks, all 
three transitions are coincident and there is no ambiguity. 
Thus, we choose to match the deconfinement transition 
at Tc = 154 MeV [33,, M], which leads us to set Tq = 50 
MeV. 

Minimizing ft in the presence of the Polyakov loop 
yields the phase diagram shown in Fig. [51 Comparing 
to Fig. [3] we now assess the effects of confinement on 
the phase structure of QCD. As in the nonconfining NJL 
model, the topology of the phase diagram depends criti- 
cally on K. We find that this dependence is unaffected by 
the inclusion of confinement, and that the critical point 
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FIG. 7: (color online) Contour plot for the chiral condensate in 
the PNJL model. 



FIG. 9: (color online) Contour plot for the diquark condensate 
in the PNJL model. 
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FIG. 8: (color online) Contour plot for the traced Polyakov loop 
in the PNJL model. 



continues to appear for 4.2 < k < 4.7, while it vanishes 
into the NG-NOR phase boundary for k > 4.8. The lo- 
cation in the phase diagram at which the critical point 
vanishes (for k = 4.8) is at marginally higher temper- 
ature in the PNJL model (T = 54 MeV) than in the 
NJL model (T = 50 MeV). This reflects the fact that 
the Polyakov loop causes a slight shift (between 2 and 4 
MeV) of the NG-NOR and COE-NOR phase boundaries 
to higher temperatures, at intermediate to high fi (note 
that the critical temperature at /i = is unchanged). 

Given that deconfinement is a high T effect, it might 
be unsurprising that it does not materially affect the low 
T critical point. However, it is important to note that 
were possible fi dependence included in the Polyakov loop 
potential, Eq. ([^5]). the present inclusion of confinement 
could have a greater effect on the phase structure of QCD. 
Unfortunately, because lattice calculations are restricted 
to fi — 0, we are unable to discern any fi dependence of 
W($,¥). 

We also note that while prior work has demonstrated 
that inclusion of the Polyakov loop pulls the NG-NOR 
and CSC-NOR phase transitions, as well as the low-T 
critical point, to higher temperatures [U, |4ll, IH] , our re- 
sults do not demonstrate such a shift. This apparent dis- 
parity is a result of the aforementioned fitting of Tq — 50 



MeV, which we choose to reproduce the /i = decon- 
finement transition for massless three-flavor QCD. Had 
we instead chosen to match the transition in the pure- 
gauge sector (as in ^2l]) and set Tq — 270 MeV, both the 
NG-NOR and CSC-NOR transitions, as well as the low- 
T critical point, would have been shifted to significantly 
higher temperatures. 

Comparing the contour plot of the chiral condensate 
(Fig.[71) to that from the NJL model (Fig.gl), we note that 
the Polyakov loop encourages larger values of tr, most no- 
tably near the phase boundaries. This is clear from the 
absence of a cr = O.Gao contour at low ^ in the PNJL 
model, as well as the termination of the a = O.Sao con- 
tour into the NG-NOR phase boundary at lower chemical 
potential (/i = 250 MeV) than in the NJL model (/i = 266 
MeV). This effect can be traced to an effective cr^$ cou- 
pling, which favors the coexistence of a chiral condensate 
and confinement, and which is discussed in more detail 
in Sec. I VII In the same vein, from Fig. [5] we note that 
$ tends to decrease in the presence of cr. As a result, 
curves of constant $ are "pulled" to higher temperatures 
in the NG phase than they would be in the absence of the 
effective (t^$ coupling. Finally, the effective quark mass 
obtained at "normal" conditions (i.e., = 0) is identi- 
cal (M — 355 MeV) whether computed with or without 
confinement. 

Inspecting the contours of the diquark condensate 
(Fig. El), we find that the Polyakov loop also has no ap- 
preciable effect on d. 

Finally, Fig. |8] demonstrates that $ is only weakly de- 
pendent on chemical potential, being primarily an in- 
creasing function of temperature, and what fi dependence 
does exist is almost entirely restricted to the NG phase. 
In the COE phase, the $ contours line up roughly par- 
allel to the d contours, suggesting an effective cou- 
pling. However, it does not appear that the coupling 
affects the deconfinement transition significantly, since $ 
has already achieved nearly its maximum value at much 
lower temperatures than where d — >■ 0. We find that 
$ is discontinuous across the first-order NG-NOR and 
NG-COE transitions (e.g. note the jump in the $ = 
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FIG. 10: (color online) Plot of the normalized chiral and di- 
quark condensates and Polyakov loop for fi = 280 McV. 




FIG. 11: (color online) Coefficient of a^, a{^,T), for various 
values of the chemical potential (solid = 0, dash = 150 MeV, dot- 
dash = 300 MeV, long-dash = 450 MeV). For a < the term 
favors spontaneous formation of a chiral condensate. 



contour at /i = 267 MeV), but the relative magnitude of 
the discontinuity is much less than that of a fFig. [TU)) . 



VI. GINZBURG-LANDAU COEFFICIENTS 

Having constructed the PNJL phase diagram, we 
next seek to understand the first-order effects of the 
condensate-Polyakov loop couplings by expanding the 
thermodynamic potential f2 in a Ginzburg-Landau (GL) 
form 



— a — - a 
2 3 



(44) 



While prior work by Hatsuda et al. focused on the topo- 
logical consequences of a thermodynamic potential of this 
form (without the Polyakov loop), here we are in a po- 
sition to compute the coefficients explicitly as functions 
of temperature and chemical potential For example, 
the coefficient a can be computed from Eq. p7p via the 
relation 



(45) 



(T=d=*=0 



while similar expressions hold for the other coefficients. 

In the following calculations, for the sake of conve- 
nience the coefficients a, c, 7, etc. will be taken to refer 
to their dimensionless versions, where they are scaled by 
the appropriate power of A (as in Table H]) . This ensures 
that the coefficients are of roughly the same order of mag- 
nitude and facilitates comparison of their relative impor- 
tance. For the purposes of these comparisons, note that 
the dimensionless chiral and diquark condensates have 
maximum values of ao = 0.0636 and do = 0.0548. In 



addition, in order to express the coefficients compactly, 
we define the following quantities: 



9^ 



1 



1 



cosh(/3i?| 
1 



1' 



(46) 
(47) 
(48) 



where = k zL fi are the eigenvalues (without absolute 
values) in the absence of any interactions. 



A. Noncoupling terms 

In computing the coefficients of simple powers of cr, 
we may set c? = $ = prior to taking the necessary 
derivatives. Thus, the 18 distinct eigenvalues shown in 
Eqs. ((38|) and ((39)) reduce to the six distinct values: 
i?i_4 = \Ek ± mI i 27rT/3 (with the two ± independent) 
and E^,Q = \Ek±^i\■ In this way, we obtain the first three 
LG coefficients: 



a(/x,r) = 12G 



(27r)3 k 
c{fi,T) = UK 



(49) 

where denotes summation over E+ and E-. We 
note that the term is proportional to the coupling K, 
which stems from the axial anomaly. This proves true 
for all odd powers of a so that in the absence of the 
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FIG. 12: (color online) (a) Coefficient of a{ii,T), for various values of the chemical potential (same as Fig. Illl l. For o < the |dp 
term favors spontaneous formation of a diquark condensate, (b) Coefficient of <&, A(/i,T), for various values of the chemical potential. 
For A{T) < the linear term favors the formation of a non-zero <&. 



axial anomaly the thermodynamic potential is an even 
function of a. 

Figure [TT] shows that the coefficient a changes sign, 
becoming negative at low temperatures and chemical po- 
tentials. As a result, in the low-/i, low-T portion of the 
phase diagram, the term in the thermodynamic po- 
tential tends to favor chiral condensation. In fact, we 
see that if the term were dominant, this transition 
would occur at /i = at the extremely high tempera- 
ture of T = 395 MeV. However, by looking at Eq. p9)) 
we can assess the relative magnitude of c, and whether 
it will play a significant role in determining the order of 
the NG-NOR phase transition. In fact, noting that the 
integrals appearing in c are identical to those in a, we 
can write 



'^=TT77(a-4G) (50) 



Thus, we find that at /i = 0, \c/a\ > 6K ~ 60, while 
as noted above, uo ^ 1/16. As expected from the results 
of the prior section then, we find that the term can- 
not be ignored in determining the order of the NG-NOR 
transition. In fact, the inclusion of higher-order terms 
coupling (T and $ leads to a first-order transition at a 
more modest temperature of Tc — 154 MeV. 

Next, the coefficients of terms involving only powers of 
\d\ may be obtained by setting cr = <I> = at the outset 
and taking the appropriate derivatives. Unlike the prior 
calculation, in which setting d = reduced the number 
of distinct eigenvalues from 18 to 6, in this case, no such 
simplification occurs and the full 18 eigenvalues must be 
evaluated. Doing so yields the |dp coefficient: 



( RU , fr2 V- / rf'k r 4sinh(^^^ 

-K9/3[6/3£;^ sinh(^£;^) + \/37r]/±/i± 
+mi3'^E^ sinh(/3i;±)5±/i±|. 

(51) 

As shown in Fig. I12al for sufficiently large /x, the sign 
of the term becomes negative at low T. Since no odd 
powers of the diquark condensate appear in i7 the phase 
transition is of second order. This transition produces a 
CSC-like COE phase of quark Cooper pairs in the low T, 
high /i portion of the pha se diagram, as has been widely 
reported [H [13, [HI, Hi ■ The precise location of the 
diquark condensate formation is affected by the chiral- 
diquark condensate coupling, but for T ~ 0, we find that 
a becomes negative at /x 120 MeV. 

Next, we move on to compute the coefficients of powers 
of $. Setting a = d = Q again reduces the 18 eigenvalues 
to six eigenvalues analogous to those involved in the cal- 
culation of the cr" coefficients (with — >■ fc). Computing 
the necesary derivatives yields: 

^(-^)-^^?/(0^^' 

S(M,r) = -r^[a(T) + 126(T)] + ^E / (^(/^)'' 
C(;.,T)=245(T)-^X:/|^(/")^- 

(52) 
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FIG. 13: (color online) (a) The coefficient of cr\d\'^, 7(ai, T), for various values of the chemical potential (same as Fig. Illl l. For 7 < 
the system favors coexistence of chiral and diquark condensates, (b) Coefficient of cr'^'J?, a'{fj,, T), for various values of the chemical 
potential. For a' < 0, the (t^<1? term favors simultaneous formation of a chiral condensate and dcconfincment. 



Significantly, while in the pure-gauge sector the lowest- 
order term in Z^($,<i>) is quadratic (see Eq. (f28| ). the 
existence of quarks generates a term linear in $. Further, 
A is negative at all points in the phase diagram except at 
T = 0. As a result, even for very small non-zero temper- 
atures, the Polyakov loop will take on a finite value. This 
behavior stands in marked contrast to that of the pure- 
gauge sector, in which $ undergoes a large discontinuous 
jump at To = 270 MeV, below which $ = 0. 



sider the lowest-order couplings between these variables. 
Beginning with the chiral and diquark condensates, we 
note from Figs. [7] and [H] that there is only a small region 
near fi ^ 280 — 320 MeV in which both a and d are sig- 
nificant. Thus, the primary condensate coupling is the 
lowest-order coupling, of the form crldp. Setting $ = 
at the outset and performing the necessary derivatives 
yields 



B. Coupling terms 



Having computed the coefficients of the pure conden- 
sate and Polyakov loop terms in Eq. (|44l) . we now con- 



7(M,T) = -i^ -2Gi^ ^1 (2/ +.g ) + ^ 

+?,HI3[I3E^ sinh(^£;^) + \/37r]/±/i± + ISHfi'^E^ sm\\{l3 E^)h^ . (53) 



In Fig. I13al we see that 7 is negative at all points in 
the phase diagram, and therefore the (7\d\'^ coupling uni- 
versally encourages coexistence of the chiral and diquark 
condensates. In the prior section we observed that this 
term is the critical factor in determining the nature of the 
NG-COE transition at low temperatures. Since the <j\d\'^ 
coupling does not involve the Polyakov loop, Eq. ([53)) 
is the same as that obtained by Abuki et al, although 



they did not compute it explicitly, but only observed its 
consequences in the numerical construction of the QCD 
phase diagram [ll|. 

Having computed the lowest-order noncoupling terms, 
we consider the effective modifications to these terms 
that arise from the inclusion of the Polyakov loop. The 
coefficients of the lowest-order couplings between the con- 
densates and the Polyakov loop, (T^$, a^^, and are: 
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(2^ k ' 



c' = 216GKJ2 I 

+3\/3/3 [2y3/3£;^ sinh(/3i;^) + vrj [f^fh^ 



(27r)3 



(54) 



As shown in Fig. ll3bl except for at very large chemical 
potential (/x > 445 MeV) both a' and c' = {iK' /2G)a' 
are positive, and therefore tend to disfavor simultaneous 
chiral condensation and deconfinement. Since a is only 
appreciable for /x < 300 MeV, this means that to lowest- 
order, the presence of a chiral condensate will tend to 
maintain a confined state, and vice versa. We note, how- 
ever, that this finding does not preclude the realization 
of a spatially inhomogeneous "quarkyonic" phase at large 
/i, as has been suggested recently p4l - l48| . since we have 
explicitly assumed a homogeneous condensate. 

We also note that the magnitude of a' decreases with 
increasing chemical potential, and is therefore largest 
(most strongly opposing deconfinement) in the NG re- 
gion of the phase diagram where a dominates. This effect 
is visible in Fig. (jS]), where the curves of constant $ are 
pulled to higher temperatures in the NG phase by virtue 
of the presence of a non-zero a. The NG-NOR transi- 
tion is therefore the result of two competing effects. On 

I 



I 

one hand, in the pure-gauge sector confinement tends to 
become weaker at higher temperatures, eventually giving 
way to a deconfined QGP. On the other hand, as tem- 
perature increases, the system has an increasing aversion 
to a state in which both a and $ are non-zero, so the 
presence of the chiral condensate tends to suppress the 
deconfinement transition. 

Similarly, as shown in Fig. I14al the coefficient a' is 
positive throughout the phase diagram and its magni- 
tude increases with increasing /i. Thus, the presence of 
a diquark condensate also tends to maintain a confined 
state with $ ~ 0, and does so more strongly at high 
chemical potentials, where d is appreciable. This can be 
understood by noting that in order to satisfy the gap 
equation, d^/dA = 0, with increasing one must de- 
crease the magnitude of the gap. 

The final coefhcient that we compute is of the term 
(t|o?P$, the lowest-order term coupling all three variables: 



/ = 6^'E/(0^{-18V3.i7/3X^ 

-A/37riJ^/±/i± 



±(/^±)2sinh(/3£;o±)5± 

SVSH^ r 



3G 



H 



sinh(/?£;±)(/±) 



±\2 



(55) 



As shown in Fig. Il4b| 7' is positive throughout the phase 
diagram, so that once again, we find that condensation 
and deconfinement tend to disfavor one another. Thus, 
we find that in general, the existence of a condensate 
(either a or d) tends to encourage confinement, and the 
greater the number or magnitude of the condensate(s) 
present, the greater the effect. 



C. Low temperature critical point 

Having computed the lowest-order Landau- Ginzburg 
coefficients, we are now in a position to assess the ef- 
fect of the Polyakov loop on the low temperature critical 
point. As noted in the prior section, in order for the 



I 

critical point to appear the ratio of the axial anomaly 
couplings K — K' /K must exceed 4.2. In the presence of 
the Polyakov loop, the values of these couplings are ef- 
fectively modified by terms proportional to $, as well as 
higher-order terms which can be neglected to first-order. 
Comparing Eqs. pT)) and pi)) we define 



Keff ^Ko 1 



(56) 
(57) 



Taking the ratio and expanding to linear order in $ yields 



1+ ( 1 $ 



(58) 
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FIG. 14: (color online) (a) Dimensionless a'{^,T) for various values of the chemical potential (same as Fig. [ni l. For a' < 0, the [dp* 
term favors the coexistence of a diquark condensate and deconfinement. (b) Coefficient of (T|iip<l?, 7'(/i, T), for various values of the 
chemical potential. For 7' < 0, the system favors the coexistence of a, d, and ^. 



The effect of ttie Polyakov loop can be assessed in terms 
of tfie sign of S = 7'/7 ~ c'/c. For 6 > 0, the effec- 
tive ratio Keff is increased, wfiicti will tend to encour- 
age the emergence of the critical point, while for 5 < 
the critical point will tend to be suppressed. The con- 
dition for Eq. ([55)1 to remain valid is (5 ^ ^ 2, 
which follows from the fact that $ ~ 0.4 at the critical 
point. Noting that for = 4.2 the critical point is at 
i^i,T) = (288 MeV, 36 MeV), we find that S^p = -0.16. 
Thus, the presence of the Polyakov loop decreases the 
effective coupling ratio n by approximately 6%. 

There are several remaining questions relevant to im- 
proving our current picture of the QCD phase diagram. 
First, the effects of realistic bare quark masses, partic- 
ularly the intermediate strange quark mass, on the low 
temperature critical point have yet to be studied in the 
context of the PNJL model. By addressing this prob- 
lem, we may come to better understand the "freezing- 
out" of the strange quark as the system moves to lower 



chemical potentials, where the light up and down quarks 
dominate. Second, we have not imposed the restrictions 
of charge-neutrality and /3-equilibrium on our analysis. 
Finally, our model might be extended by including ad- 
ditional condensate structure beyond CFL, which would 
allow for a diversity of CSC-COE phases in the low tem- 
perature portion of the phase diagram. 



VII. ACKNOWLEDGEMENTS 

This research was supported in part by NSF grant 
PHY09-69790. The authors thank Professor Tetsuo Hat- 
suda for insightful comments throughout the course of 
this work, and for critically reading the manuscript. P.P. 
also thanks Tomoki Ozawa for posing helpful questions 
and perspective. 



[1] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 
(1961). 

[2] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 124, 246 
(1961). 

[3] M.K. Volkov and A.E. Radzhabov, Phys. Usp. 49, 551 
(2006). 

[4] R.D. Pisarski and F. Wilczek, Phys. Rev. D 29, 338 
(1984). 

[5] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 44, 
1422 (1970). 

[6] G. 't Hooft, Phys. Rev. Lett. 37, 8 (1976); Phys. Rev. D 

14, 3432 (1976). 
[7] D.J. Gross, R.D. Pisarski, and L.G. Yaffe, Rev. Mod. 

Phys. 53, 43 (1981) 



[8] J.K. Boomsma and D. Boer, Phys. Rev. D 81, 074005 
(2010). 

[9] T. Hatsuda, M. Tachibana, N. Yamamoto, and G. Baym, 

Phys. Rev. Lett. 97, 122001 (2006). 
[10] N. Yamamoto, M. Tachibana, T. Hatsuda, and G. Baym, 

Phys. Rev. D 76, 074001 (2007). 
[11] H. Abuki, G. Baym, T. Hatsuda, and N. Yamamoto, 

Phys. Rev. D 81, 125010 (2010). 
[12] T. Hatsuda and T. Kunihiro, Phys. Rep. 247, 221 (1994). 
[13] M. Buballa, Phys. Rep. 407, 205 (2005). 
[14] T. Schafer, Proc. Hampton Univ. Grad. Stud. Prog, (to 

be published). 

[15] Chr. V. Christ ov, K. Goeke, and M. Polyakov, 
|hep-ph/950i383l (1995). 



13 



[16] H.J. Warringa, D. Boer, and J.O. Andersen, Phys. Rev. 

D 72, 014015 (2005). 
[17] J.K. Boomsma and D. Boer, Phys. Rev. D 80, 034019 

(2009). 

[18] K. Fukushima, Phys. Lett. B 591, 277 (2004). 
[19] C. Ratti, M.A. Thaler, and W. Weise, Phys. Rev. D 73, 
014019 (2006). 

[20] C. Sasaki, B. Friman, and K. Redhch, Phys. Rev. D 75, 
074013 (2007). 

[21] S. Rofiner, C. Ratti, and W. Weise, Phys. Rev. D 75, 
034007 (2007). 

[22] P. N. Meisinger and M. C. Ogilvie, Phys. Lett. B 379, 
163 (1996). 

[23] H. Abuki, Prog. Theo. Phys. Supp. 174, 66 (2008). 

[24] A. M. Polyakov, Phys. Lett. 72B 477 (1978). 

[25] E. Megi'as, E. R. Arriola, and L.L. Salcedo, Phys. Rev. 

D 74, 065005 (2006). 
[26] T. Banks and A. Ukawa, Nucl. Phys. B 225, 145 (1983). 
[27] J. Kogut et al, Phys. Rev. Lett 50, 393 (1983). 
[28] T. Z. Nakano, K. Miura, and A. Ohnishi, Phys. Rev. D 

83 016014 (2011). 
[29] A. Bazavov and P. Petreczky, J. of Phys. Conf. Ser. 230 

012014 (2010). 

[30] T. Schafer and F. Wilczek, Phys. Rev. Lett. 82, 3956 
(1999). 

[31] B. Svetitsky, Phys. Rep. 132, 1 (1986). 
[32] L.G. Yaffe and B. Svetitsky, Phys. Rev. D 26, 963 (1982). 
[33] K. Fukushima and T. Hatsuda, Rep. Prog. Phys. 74, 
014001 (2011). 

[34] F. Karsch et al., Proc. Strong and Electroweak Matter 

Meeting (to be pubhshed). 
[35] Y. Aoki et al., Phys. Lett. B 643, 46 (2006). 
[36] Y. Aoki et al., J. of High Energy Phys. 6, 88 (2009). 
[37] There is a shght ambiguity in this definition at low /i. 

As shown in Fig. 1101 for fi < 288 MeV, as the system 



moves to higher T, $ first undergoes a smooth, though 
rapid, crossover, which is followed by a slight discontinu- 
ity coincident with the chiral phase transition. Whether 
one chooses to identify the deconfinement transition with 
the maximum finite value of d^/dT (rapid crossover) or 
with the absolute maximum of d^/dT (discontinuity) is a 
matter of convention. We choose the latter definition. For 
large fi there is no such ambiguity, since the COE-NOR 
transition is second order and no discontinuity exists. 
[38] F. Marhauser and J.M. Pawlowski, hep-ph/0812.1144 
(2008). 

[39] J. I. Kapustaand C. Gale, Finite- Temperature Field The- 
ory, 2nd ed.. New York: Cambridge University Press, 67 

(2006) . 

[40] T. Gjestland, Gauge theories at finite temperature (Mas- 
ter's Thesis), Norwegian Univ. of Science and Tech. 

(2007) . 

[41] D.B. Blaschke, F. Sandin, V.V. Skokov, and S. Typel, 
Acta Phys. Pol. B Proc. Suppl. 3, 741 (2010). 

[42] D. Gomez Dumm, D.B. Blaschke, A.G. Grunfeld, and 
N.N. Scoccola, Phys. Rev. D 78, 114021 (2008). 

[43] M. Alford, Ann. Rev. Nucl. Part. Sci. 51, 131 (2001). 

[44] L. McLerran and R. D. Pisarski, Nucl. Phys. A 796, 83 
(2007). 

[45] Y. Hidaka, L. D. McLerran, and R. D. Pisarski, Nucl. 

Phys. A 808, 117 (2008). 
[46] K. Fukushima and T. Hatsuda, Rept. Prog. Phys. 74, 

014001 (2011). 

[47] S. Carignano, D. Nickel, and M. Buballa, Phys. Rev. D 

82, 054009 (2010). 
[48] T. Kojo et al, Nucl. Phys. A 875, 94 (2012). 



